Astronomy & Astrophysics manuscript no. AA_2009_1 3649v5 


©ESO 2010 


November 15, 2010 





Black Holes and Galactic Density Cusps 
From Black Hole to Bulge 

M. Le Delliou 1 , R.N. Henriksen 2 , and J.D. MacMillan 3 



o 

> 
O 

(N 



Instituto de Ffsica Teorica UAM/CSIC, Facultad de Ciencias, C-XI, Universidad Autonoma de Madrid 
Cantoblanco, 28049 Madrid SPAIN 
e-mail: Morgan.LeDelliou@uam.es 

2 Queen's University, Kingston, Ontario, Canada 
e-mail: henriksnOastro . queensu . ca 

3 Faculty of Science, University of Ontario Institute of Technology, Oshawa, Ontario, Canada L1H 7K4 
e-mail: j oseph . macmillan@gmail . com 

Abstract 

Aims. In this paper we continue our study of density cusps that may contain central black holes. 

Methods. We recall our attempts to use distribution functions with a memory of self-similar relaxation, but mostly they apply only in 
restricted regions of the global system. We are forced to consider related distribution functions that are steady but not self-similar. 
Results. One remarkably simple distribution function that has a filled loss cone describes a bulge that transits from a near black hole 
domain to an outer 'zero flux' regime where p <x r~ 7//4 . The transition passes from an initial inverse square profile through a region 
having a l/r density profile. The structure is likely to be developed at an early stage in the growth of a galaxy. A central black hole is 
shown to grow exponentially in this background with an e-folding time of a few million years. 

Conclusions. We derive our results from first principles, using only the angular momentum integral in spherical symmetry. The initial 
relaxation probably requires bar instabilities and clump-clump interactions. 

Key words, theory-dark matter-galaxies:haloes-galaxies:nuclei-black hole physics-gravitation. 
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We discussed the relation between the formation of black holes 
(hereafter BH) and of galactic bulges, or of bulges with a central 
mass in the previous papers of this series of papers (Henriksen 
et al., referred to here as paper|T]|II] and references hereafter). 

Previous research (IKR19951 IMa98l) and more recently 
■ dFerrase and Merritt 20001 IGebhardt et~al"20 00 ) has established 
a strong correlation between the BH mass and the surround- 
ing stellar bulge mass (or velocity dispersion), which we take 
as an indication of coeval growth. Such growth may occur as 
BH 'seeds' accrete gas or disrupted stars in a dissipative fash- 
ion during the AGN (Active Galactic Nuclei) phase, but there 
is as yet no generally accepted scenario. Moreover, the ne- 
cessity for a seed and observations of early very supermas- 
sive BHs (e.g. |Kurk et al., 2007| l, together with changes in the 
normalization of the BH mass-bulge mass proportionality (e.g. 
Mai olino et al. 2007b : suggest an alternate growth mechanism. 
In fact some authors dPeirani & de Freitas Pacheo 20 08) have 
studied the possible size of the dark matter component in BH 
masses and deduce that between 1% and 10% of the black hole 
mass could be due to dark matter. 

We explored coeval growth in the case of spherically sym- 
metric radial and non-radial infall (paper [I] [II]), using distribu- 
tion functions with a self-similar memory. However although we 
obtained reasonable descriptions in sections, we did not find a 
distribution function (DF) that could be fit to Black Hole halo 
and bulge continuously. Moreover the problem of the black hole 
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feeding was left unresolved. This requires an anisotropic DF 
with a filled loss cone of the nature of the Bahcall & Wolf solu- 
tion dBahcall & Wolf 1976b . 

Here we repair these omissions for a special case of non 
self-similar infall. We use the same technique of inferring rea- 
sonable distribution functions for collisionless matter from the 
time dependent Collisionless Boltzmann (CBE) and Poisson set, 
while insisting on a central point mass. The temporal evolution 
allows for relaxation of collisionless matter in addition to possi- 
ble 'clump-clump' (two clump) interactions. 

We use, in this paper as in the previous (|I] |nj, the Carter- 
Henriksen dCarter & Henriksen 1991b procedure. In this way 
we obtain a quasi-self-similar system of coordinates (H2006 
H2006A) that enables expression of the CBE-Poisson set with 
explicit reference to a previous transient self-similar dynamical 
relaxation. Thus we can remain 'close' to self-similarity just as 
the simulations appear to do. 

We begin the next section with a summary of this approach 
in spherical symmetry that is common to all papers in this se- 
ries, including the useful results found earlier (papers [I] [TIJ. 
Subsequently we recall the limiting case with a high binding en- 
ergy cut off in the distribution function (hereafter DF). This is to 
remind us that such a cut-off, existing for any reason, produces 
a flat cusp. Finally we give the key result of this paper as a non- 
self-similar DF (of a type found previously dEvans & An 2 006 ) 
but not used in this connection) that describes the transition zone. 
This zone extends from a central mass through a shallow re- 
gion (to be compared with observations of the milky way) into 
a Bahcall & Wolf feeding region. Ultimately it enters an inverse 
square law region. 
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2. Dynamical Equations and radial results 

Following the formulation of H2006 we transform to infall vari- 
ables the collisionless Boltzmann and Poisson equations for a 
spherically symmetric anisotropic system in the 'Fujiwara' form 
(e.g. Fujiwara 1983) namely 



dt 
d 



df 
or 
2 d<P 



r 3 dr J dv r 



— A% G 



j f(r,Vr,j 2 )dv r df 



(1) 

(2) 



where / is the phase-space mass density, <t> is the 'mean' 
field gravitational potential, j 2 is the square of the specific angu- 
lar momentum and other notation is more or less standard. 

The transformation to infall variables has the form (e.g. 
IH2006I I 



aT 



R = re- aT l a , 

Z = / e -(4/«-2)«r e -* =at 

P(R,Y,Z;T)=e^ a - 1 '> aT nf(r,v n j 2 ;t), (3) 
y{R;T)=e- 2{l /"- l)aT <i>(r), 
®(R;T)=p(r,t)e- 2aT . 

The passage to the self-similar limit requires taking dj = 
when acting on the transformed variables. Thus the self- 
similar limit is a stationary system in these variables, which is 
a state that we refer to as 'self-similar virialisation' flHW 1999 
ILe De lliou 2001). The virial ratio 2AT/jW| is a constant in this 
state (although greater than one; K is kinetic energy and W is 
potential), but the system is not steady in physical variables as 
infall continues. 

The single quantity a is the constant that determines the dy- 
namical similarity, called the self-similar index. It is composed 
of two separate reciprocal scalings, a in time and 8 in space, 
in the form a = a/ 5. As it varies it contains all dominant phys- 
ical constants of mass, length and time dimensions, since the 
mass scaling jj. has been reduced to 35 — 2a in order to main- 
tain Newton's constant G invariant (e.g. H2006). 

We assume that time, radius, velocity and density are mea- 
sured in fiducial units r /v ,r , v and p„ respectively. The unit 
of the distribution function (DF from now on) is f and that of 
the potential is v 2 . We remove constants from the transformed 
equations by taking 



fo = Po/vl, 



\%Gp r L . 



(4) 



These transformations convert equations to the re- 

spective forms 



id T P-(3/a-l)P+(---)d R P 
a a a 

(l/ a -l)Y+U^-^\\d Y P-(4/a-2)Zd z P = 



and 



R 2 dR \ dR 



= ©. 



(5) 
(6) 



This integro-differential system is closed by 
(-) I PdY <IZ. 



(7) 



This completes the formalism that we will use to obtain our 
results. 

For easy reference in this paper, we summarize the general 
forms of the self-similar DFs that we discussed previously. 

In paper [T] we discussed the emergence in purely radial in- 
fall of the rigorously steady DFs from Henriksen and Widrow 
(hereafter HWDF, IHW95b 



nf = F(K)\E\ l / 2 SU 2 ), 



(8) 



and the time-dependent DF from Fridmann and Polyachenko 
(hereafter FPDF, Fridman n~& Polyachenko 1984] ) 



K 



-E+E yi 2 



8(f)- 



(9) 



Introducing anisotropics in paper [II] we found the general self- 
similar solution to equation (|5]l to be 



P = P(f,K)£«, 



where 

l -a 

K = SZ'tt 
3 — a 
q = 2(o-l) 



(10) 

(ID 
(12) 



This yields the physical form of the general steady state with 
self-similar memory to be 



K 



P(K)\E\«, 
\E\(f)-l& 



(13) 



where E = v 2 /2 + j 2 / (2r 2 ) + <f>. Taking first P to be constant, 
and then taking it to be proportional to K~ q . this yields respec- 
tively 



nf 
nf 



K\E\i, 
K 



(14) 
(15) 



where we have defined w = (3 — a)/(4 — 2d). These limits were 
found to apply to inner and outer extremes of the relaxed region 
of a simulated 'bulge' (i.e. no black hole was included in the 
simulations). 

The DF ( TT3l can be put in a form that seems to generalize 
the FPDF (Eq.|9]). We choose P« rc~(^) to obtain 



nf-. 



K 



(j 2 )^\E -E\V 2 



(16) 



The density and potential laws are more general than those of 
the FPDF, being respectively p °c r~ 2 " and <t> oc r 2{ - x ~ a \ 

None of these DFs can be said to apply to every region of the 
simulated halos. 
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3. High Binding Energy Cut-off 

We recall in this section an extreme case that limits the 
flattening of the cusp near a black hole. There is evidence 
dLe Delliou 20011 IMWH 20061 1 for a cut-off in the number of 
particles with high binding energy, even without a central point 
mass. The presence of a central black hole of mass M. is likely to 
accentuate this trend by accretion. We may imitate such a cut-off 
in order to study limiting behaviour near the black hole, by using 
the isothermal DF (see e.g. papers I and II) at negative energies 
and temperatures in the form 



nf = Ke^ E \ 



(17) 



Here E = v 2 r /2 + j 2 /(2r 2 ) +4> < 0. We suppose that the con- 
stant j3 > and corresponds to some reciprocal cut-off energy. 
A straight-forward calculation of the density implied by such a 
DF yields (note that we integrate only over negative energies 
rather than over all velocities) 



8V2K 



<I>| 3/2 M(l,5/2,-/3|4>|), 



(18) 



Where M(a,b,z) is the Kummer function. We expect this to 
hold where \<S>\ sw GM,/r is large, and in this limit the Kummer 
function varies as — 3/(2z). Consequently the density cusp be- 
comes as flat as r~ l l 2 . This was already noticed in 

dNakano & Makino 1999t and in (IMerritt & Szell 20061 1. where 
in both cases it is due to scouring by merging black holes. In any 
case, as was recognized in (Nakano & Makino 1999), only the 
cut-off is required. 

A cut-off in the number of particles per unit energy dM/dE 
is to be expected on general grounds. The density of states 
(Binney & Tremaine 1987 ) for an isotropic DF in a region dom- 
inated by a point mass potential is g(E) — V2(7tGM,) 3 \E\~ 5 / 2 , 
from which the differential energy distribution dM/dE = 
g(E)f(E) may be calculated. Thus even a constant f(E) (for 
which p oc |<J>| 3 / 2 and therefore goes like 5 near a point mass) 
will show a cutoff in the mass distribution as the density of states 
decreases in the expanding phase space. However the evidence 
suggests that f(E) is also declining with increasing \E\ in the 
most tightly bound central regions. 

Assuming the DF (TBI) and a — 0.5), we predict the mass 
distribution function dM/dE °< |£| -5 . Thus a cut-off in the 



isotropic DF can be associated with the r = r density 
profile. The mechanism for the high binding energy cut-off is 
not immediately evident, but it must be part of the dynamical 
relaxation. High negative energy particles must be preferentially 
excited to less negative energies, which does not seem to oc- 
cur in strict shell code simulations dHW 19991 1 . However it may 
occur due to the presence of sub-structure such as clumps or 
bar formation by the radial orbit instability (MWH 2006J. In 
any case it may not require black hole 'scouring', although this 
process does have the maximum effect (Merr itt & Szell 2 006, 
INakano & Makino 1999ft . 



4. Global Solution with Black hole and a 
Bahcall/Wolf Outer Cusp 

We consider an exact collisionless cusp that has an embedded 
central mass. It is necessarily not self-similar. We have seen that 
we expect a DF that is primarily dependent on angular momen- 
tum in the outer regions of a bulge (e.g. ([TBI ) before the maxi- 
mum in angular momentum). We can not follow analytically the 



development of this region in time, as this is the province of nu- 
merical work. We can however hope to find analytic descriptions 
of ultimate states. 

There is one steady solution that is not self-similar but which 
is closely related to the Evans and An (Eva ns & An 2 006) solu- 
tion with anisotropy parameter j3 = 0.5. This is an exact solution 
to equation (0 with dr ^ (compare equations Q7J [TBI ) in the 
form (recall that w = (3 — a) /{A — 2a)) 

p = p(f)z- n \ 

that is, written entirely in transformed coordinates, 



P = P(Zexp (A/a - 2)aT)Z~ 



(19) 



This solution is always steady but not always self-similar be- 
cause of the T dependence. This equation is seen to satisfy equa- 
tion © by direct substitution when the dependence on T is in- 
cluded. The steady form ( TTOb that we found in the previous paper 
of this series is distinguished from the general steady class by be- 
ing a possible limit of self-similar evolution. The same problem 
has been addressed with arbitrary steady and isotropic distribu- 
tion functions in an interesting paper by Baes et al. (2005) 

A convenient choice is to take P — K(j 2 ) v , where v is any 
suitable real number. Then the DF is simply 



nf = K[f) 



2\(v-w) 



(20) 



The density associated with such a distribution is (E < 0, 

V— w> — 1) 



r(v-i 



n 



(v- 



■3/2)r(v- 



w - 



•3/2) 



„2(v— vv) | 2<J> |(V— vv+3/2)_ 
(21) 

This may now be substituted into the Poisson equation to give an 
equation for |<f>| and hence the density in a self-consistent bulge. 

There is a power law solution °c r h provided that h > — 1 . In 
terms of w and v, h is 



h = 



-2(1 



• w) 



v-w + 1/2 ' 



(22) 



and the logarithmic density slope is h — 2. Most of the density 
power laws (h — 2) of interest arise for values of v — w near — 1 . 
For example v — w — —0.9 yields h — 2 = —1.5, while 

V — w = —5/6 yields h — 2 = —1.0. The limiting case with 

V — w — — 1 (which has a logarithmic singularity at zero angular 
momentum) gives as expected h — 2 = —2. Setting v = returns 
us to the self-similar bulge ( PTBT ) and indeed h — 2= —2a in this 
case. 

It is clear however that there is a special case when v — w + 

1/2 = 0, for which the DF is nf = K/ \/~j 2 ~. This value treated 
separately yields 

P = — 1*|. (23) 
r 

The Poisson equation is therefore linear and it is readily solved 
to find 



1*1 = 



(24) 

where J\ and Y\ are first order Bessel functions of the first and 
second kinds respectively. The density is given by equation ( l23l . 

The constants in equation d24l) must be such as to maintain 
|<£>| and p positive despite the oscillations in the Bessel func- 
tions. This may impose a limited range in radius, which is in 




Figure 1. We show the various density regimes of the zeroflux cusp with an embedded black hole. The constants in equation 
d24t have been chosen for numerical convenience and positive definiteness to be: |4>|«, = 3, Ci = 5 and C2 = —2/(nK'). Here 
K' = K/(8k) — 10 so that C2 = —M, = — l/5n. In addition the plot is of p / \fW . All quantities are expressed in fiducial units; r a , 
po and v c . 



any case required to ensure a finite mass. In general the Bessel 
functions will be out of phase by n/2 and it is not possibile to 
keep their sum positive everywhere. However it is possible to 
keep the whole expression positive by choosing an appropriate 
set of constants as we have done in our example. Because the 
term containing the Bessel functions is declining as r~ 7//4 , once 
a critical transition region is traversed successfully by a choice 
of constants the positivity is assured henceforward. It is really 
only the difference of the potential moduli that is physical and 
this difference is free to oscillate as indicated. It is interesting 
nevertheless that |4>|oo is necessarily non-zero, as this suggests 
that the solution must be embedded in distant matter. 

As z — > one has J\ (z) » z/2 while Y\ (z) —2/{%z). Hence 
the potential contains a central point mass in this limit if C2 = 
— M„. The term in J\ tends to a constant and is presumably asso- 
ciated with the bulge mass itself. Consequently by equation ( 1231 
the density cusp near the central black hole has the profile 



p=2nK 



ci + m, 



GM. 



(25) 



The constant C\ gives the difference in the modulus of the bulge 
potential, coming from infinity to a radius where J\ is well ap- 
proximated by the small argument limit (it depends on K). This 
happens before the small argument limit applies to Y\ . We see in 
general from this solution that the near black hole cusp can start 
with a (—1) logarithmic density slope and steepen to (—2) near 
the black hole. This appears to encompass the observations. 

A remarkable behaviour of this solution occurs at large r 
(actually large y/SnKr), where Ji (z) x y/ 2 /(nz) cos (z - 3n/4) 
and Y\(z) x ^/2/(^z)sin (z— 2>k/A). This shows that the mean 



density tends to the r 7 / 4 density profile that is characteristic of 
the Bahcall and Wolf dBahcall & Wolf 19761 1 zero flux solution. 
According to the small r limits, we find it here outside a flat- 
ter region r _1 that gives way ultimately to an inner r~ 2 region. 
Outside of the Bahcall/Wolf r~ 7 / 4 region one would return to a 
r profile (recalling that the potential there, i.e. at 'infinity' is 
not zero). This appears to connect a black hole through interme- 
diate cusps to the inner NFW profile. Our figures are illustrative 
however and the constants of our formula would have to be cho- 
sen for the best fit consistent with positivity in any particular 
case. 

Figure (Q3 indicates these various regimes when the con- 
stants are chosen such that the bulge mass is about fifteen (ac- 
tually 5n, we represent the bulge mass by M a = 4np rl) times 
that of the central black hole. Starting from upper left the curves 
move inwards in a clockwise direction to show the outer r~ 7 / 4 
oscillating region, then the r region, a transition region, and 
finally the r region. The curves are connected in a continuous 
fashion. 

The final figure (|2|l expands the outer oscillating Bahcall- 
Wolf region. 

Because of the slow decline in the density with radius, it is 
evident that the outer radial limit to this solution is finite. Exactly 
where it applies will depend on special numerical cases, but it 
seems in any case that it must be inside the radius where the 
density has an inverse square profile, essentially the scale radius 
of the dark matter simulations. 

The physical behaviour all stems from a DF that has the 
simple form K/\j\. It appears to describe a zero flux equilib- 
rium condition in the mean, together with a stably filled loss 
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Given the e-folding time on the right of this last equation, we 
can increase the mass of the black hole one hundred fold in 
three to four million years. On this basis one can expect the 
presence of massive black holes very early in the history of the 
Universe. This has to happen however before the bulge has be- 
come isotropic, presumably in the rapid growth phase of the dark 
matter halo. Some seed mass is nevertheless required, which 
may be primordial. 

Such a multi-power law behaviour does not give the sim- 
ple linear correlation between bulge mass and black hole mass 
that we suggested in papers I and II. However since the black 
hole is incorporated into the global solution at all times, we do 
expect such a correlation from the integration of equation (T231 
integrated over r. If most of the inner part of the outer r~ 7 / 4 
law can be considered to have fallen into the black hole, then 
M./M b » {r./r b fl\ 



Figure 2. This figure shows the oscillations in the outer zero 
flux regime. 



cone (Tremaine 2005 The density oscillations may indicate that 
collective behaviour is necessary. It is an isolated example of 
the class of distribution functions studied by Evans and An 
( Evans & An 20061 

We can estimate the growth of a black hole in such a bulge 
in a straightforward way according to 



dM. 

dt 



= 4^p.(v r )., 



(26) 



where r. = 6GM,/c 2 . The required mean radial velocity is found 
to be ~ r — A\/2K/ (3pr) |4>| 3//2 which is the same for the inward 
or outward going particles in a steady state. We calculate v r from 
the expression 



v r = -75— I v r dv r 
r 2 p 



df 



(; 2 ) 1/2 ' 



(27) 



which becomes explicitly on transforming to an energy variable 
through dE = v r dv r 



K 



r A p J<s> 



dE 



2r 2 (£-#) 



df 



(,-2)1/2- 



This is conveniently written as 
2K 



r 2 p 



3- d\E\J2r 2 (\<P\-\E\) 



(28) 



(29) 



which yields finally the quoted result. In effect we have by ig- 
noring dE = —v r dv r calculated only the outward mean velocity 
which we take to be equal to the inward mean velocity of par- 
ticles reaching the black hole. We eliminate K in terms of the 
bulge quantities using equation ( 1231 so that 2nK = ri,Pb/\^\b- 
By setting |<t>|. = c 2 /6, the calculation yields finally (restoring 
units) 

dlnM. 8 f r b GpbC^ 



dt 



1*1, 



(30) 



As an illustration we take p b to be the mean density of 10 9 M© 
in one kiloparsec (i.e. r b ) and set |<£|j = v 2 irc (rb), the circular 
velocity at r b . We take this to be 200 km/sec and obtain finally 



iilnM. 

dt 



1.3 x lQ- 6 (y r y 



(31) 



5. Conclusions 

We have sought in this series of papers ( I II and the present pa- 
per) to find distribution functions that describe both dark matter 
bulges and a central black hole or at least a central mass con- 
centration. In most cases we succeed only in describing the dark 
matter bulges, but there are some notable exceptions. 

In the discussion of cusps and bulges based on purely ra- 
dial orbits (paperQ), we were able to distinguish the Distribution 
function of Fridmann and Polyachenko (0 from that of 
Henriksen and Widrow |(8). The FPDF was found to describe 
accurately the purely radial simulations of isolated collision- 
less halos carried out in (MacMillan 2006). These simulations 
retained the initial cosmological conditions although non-radial 
forces were switched off. The final state is close to self-similar 
virialisation rather than steady virialisation, since the infall con- 
tinues. 

This correspondence between the theory and the simulations 
gives us some confidence in the DF's found by remaining 'close' 
to self-similarity. This is especially so since predictions describ- 
ing the simulation results in paper II were based on a — 0.72 , 
which was deduced elsewhere in the context of adiabatic self- 
similarity JH2007I I. 

The FPDF can contain consistently a central mass concen- 
tration that is unlikely to be a true black hole, at least in the early 
stages. It may represent a central mass concentration or bulge 
initially. Subsequently with the rise of dissipation and instabili- 
ties, there may be a slower phase of radial accretion towards the 
centre. It is possible that this cycle could repeat several times in 
a process we have referred to as 'interrupted accretion' . Under 
this process the r~ 2 density law would apply almost everywhere. 
The radial velocity dispersion is proportional to the potential. 
Thus it decreases as near the central mass, and subsequently 
decreases logarithmically with r. 

The HWDF ([8]) is restricted to a strictly steady and self- 
similar bulge, but it has the merit of allowing a family of 
densities and potentials (velocity dispersion) according to the 
self-similar prescription. A central mass is allowed only in the 
Keplerian limit wherein a = 3/2. This gives a massless bulge 
with p r~ 3 . This is naturally iterated to give an inner flatten- 
ing but continued iteration is effectively in powers of lnr, which 
should therefore yield weak corrections. The iteration can only 
apply away from r = 0, so that the central mass is a 'renormal- 
ized' extended mass. 

The HWDF was shown in paper I to give a density that is 
linear in the potential, and hence a self-consistent bulge is found 
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from the Poisson equation. The density profile is never flatter 
than r~ 2 5 near the central mass and tends to r~ 3 in the near 
Keplerian limit of dominant central mass. This restricts the ap- 
plicability to a region outside the central bulge. It does not seem 
to be relevant to a near black hole domain. 

The inclusion of angular momentum led to more realistic sit- 
uations. We re-derived in paper ITTI the steady self-similar DFs 
from first principles in equation CfT~3T > . We showed that these can 
be used to describe the simulated collisionless halos calculated 
in (MacMillan 2006 ), if we use the value of a m 0.72 identified 
in dH2007l l and take two limits. In one dT4l > the DF is isotropic 
and describes approximately the central region of the bulge. The 
other limit (TT3T > describes the outer region. This encourages us re- 
garding the relevance of this family. Unfortunately the potential 
of a central mass can not be included exactly in these distribution 
functions. 

As always a = 1 must be treated separately and we give a 
derivation from first principles in this series. In paper I we show 
that this case corresponds to a radially growing system with the 
FPDF. The anisotropic solution is discussed in paper II The re- 
sults are new. It implies the density profile p °c r always. 

Generally we find that we can not describe elegantly 
anisotropic bulges containing black holes with self-similar DFs, 
as one might expect. The self-similarity is restricted to the sur- 
rounding bulge, and this is progressively perturbed near fhe^ 
black hole. We re-emphasized in the current paper in a non self- 
similar context that a sharp cut-off of any kind (whether due to 
black hole binary scouring or not) can yield a density profile as 
flat as -1/2. 

Our most successful description of a black hole embedded 
in an anisotropic bulge is given in this paper by the DF nf = 
K/(j 2 yl 2 . This simplest member of the family discovered by 
Evans and An (ibid) yields a bulge containing a central point 
mass. The density profile near this mass is given by equation 
( |25] | as with an inner r~ 2 peak. Farther out there is an r~ 7 / 4 
domain in the mean that imitates the Bahcall and Wolf (ibid) 
cusp, which in turn ultimately becomes the NFW r~ 1 profile. 

The reasons for this behaviour appear to be very different 
from those of Bahcall and Wolf, since there are no two body 
collisions in this treatment. By assuming a steady filled loss cone 
of this simple form, we have selected a zero flux behaviour on 
the outer boundary as an average behaviour independently of 
the detailed mechanism that acts to establish this (df/dE = 0). 
In that sense it is a useful description of the black hole cusp 
region. The oscillations may indicate the necessity of collective 
behaviour for maintaining a filled loss cone. 

We observe finally that this DF does not require a very strong 
concentration of low angular momentum particles. In the outer 
region where <t> w constant, we find dN/dj 2 ~ constant so that 
the particles are mostly at high angular momentum, as one might 
expect. 

Black hole growth is astrophysically rapid in this distribu- 
tion as seems to be required observationally. Careful study of 
the density regimes may permit bulge mass and black hole mass 
to be distinguished and correlated, assuming this distribution is 
realized. 
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